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SECTION  I 
INTRODUCTION 


A concise  mathematical  model  for  the  representation  of  the  flow  about 
two  or  mote  axisymaetric  bodies  at  zero  angle  of  attack  in  irrotational, 
incompressible  flow  has  been  presented  by  Martin  (Reference  1,2).  The  purpose 
of  this  investigation  is  to  modify  their  method  to  take  into  account  free 
stream  cross  -flow  due  to  angle  of  attack  and  to  allow  the.  consideration  of 
slightly  non-axisymmetric  bodies. 

The  method, as  given  in  Reference  2, consists  of  distributing  a set  of 
point,  sources  along  the  axis  of  each  body  which  represents  that  body  in  an 
isolated  flow  situation, and  then  adding  to  this  a set  of  source-sink  pairs 
that  preserves  the  body  shapes  when  the  bodies  are  brought  into  close  prox- 
imity. The  method  uiscussed  herein  retains  the  system  of  sources  and  source- 
sink  pairs,  but  adds  to  it  a system  of  doublets  that  accounts  for  both  the 
freestreara  cross-flow  due  to  angle  of  attack  and  the  flow  due  to  distortion 
of  the  body  centerlines. 

The  new  technique  was  applied  to  calculate  the  flow  field  about  one, 
two,  and  three  M-117  stores.  The  results  of  these  calculations  are  pre- 
sented in  Section  III. 


SECTION  II 

theoretical  considerations 


INTRODUCTION 


-r  'lnf '£'««  eVal"ati“"  °f  • 

bodies  in  close  proximity  1„  IrrotaM™  ? . ar°“nd  tuo  »r  »°re 

strict inn  is  made  that  the  bod^^s'^e^eiJhe^vr"68811516  floW‘  ^ re~ 
lution,  or  bodies  which  are  derive/f*-  * H v axisynBBetric  bodies  of  revc- 
0r  curving,  the  centerline  Also  th  °f  revol«tion  by  warping 

™«t  be  small . A1S0’  the  an8^  of  attack  of  the  bodies 


lion  2 r^ZZtT.:'  in  the  Upi.ce 


V2*  * 0 

vith  the  usual  body  boundary  condition 


(1) 


expressed  by 


V*grad(f)  - o 


(2) 


Note  that  f(x.y.z)  . n *<,  -u 

linear,  the  method  of  superposition^  °fc  ^ body*  Slnce  Equation  (1)  is 
as  given  in  Reference  2 fo^tu^a! U*ed‘  Therefore/the  solution 
axl8y “metric  bodies  can  be  e«!!S  °f  tWO  si“*lar 

of  f°  69 ’ Als°’  by  the  addition  of  doublet/*1  J°  /clude  a greater  number 
of  freestream  cross-flow  can  be  taken  into  a thelr  lna*ea’  the  effec ts 
that  the  mutual  interference  effects  oJ  tL  /!/  ^ mthoi  assuines 

""  " - 8h6d  ^ 
foiled  "UltiPle  **  “l«i-  in  generated  1„  four  b.slc  8tep8  „ 


line.  Thetr  IrlVdAtl  TT^  “•«***  -Ion,  the  c.„„ 

treating  the  body  as  an  axls^LJh?  ’’k'*!!  b°""<1*r>'  condition! 
with  a velocity  equal  to  the  axial  cornel  b°dy  allP>ed  in  a flow 

of  point  Bources  will  be  referred  to  es  the  elLl^h^"8"""-  1Ms  8et 

tne  single  body  sources. 

2.  Next,  for  each  body  a spt  nf  j v. 
with  the  single  body  sources,  represent^  t8  iB.£?md  vhich>  combined 
Includes  both  the  warped  bod^  effect!  * A a'tual  b°dy  shaPe- 

angle  of  attack.  7 eItects  "d  th*  freestream  cross-flow  due  to 
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3.  In  the  third  step,  a set  of  image  source-sink  pairs  for  each  body 
are  found  for  each  possible  combination  of  bodies  taken  a pair  at  a time. 

This  step  is  similar  to  the  method  as  given  in  Reference  2,  but  takes  into 
account  the  warp  of  the  bodies. 

4.  Finally,  sets  of  image  doublets  are  found  for  each  unique  pair  of 
bodies.  The  body  boundary  conditions  used  in  this  analysis  take  into  account 
the  flow  field  generated  by  all  sources,  image  pairs,  and  doublets  which 
were  previously  found  for  the  two  bodies  of  the  pair. 


DISCUSSION  OF  BODY  GEOMETRY 

A basic  body  of  revolution  is  geometrically  described  in  three  dimen- 
sional space  by  expressing  the  local  body  radius,  R,  as  a function  of  axial 
position,  x.  That  is,  R is  an  explicit  function  of  x where  x ie  measured 
along  the  axis  of  revolution  from  nose  to  tail  as  shown  in  Figure  1(a).  The 
non-axi3ymmetric  bodies  used  in  this  work  are  obtained  by  taking  the  body  of 
revolution  and  distorting  it  by  warping  the  axis,  as  shown  in  Figure  1(b), 
while  the  local  cross-section  taken  perpendicular  to  the  new  axis  remains 
circular  and  of  the  same  radius  and  slope  relative  to  the  centerline  as 
the  cross-section  at  the  corresponding  station  of  the  unwarped  body.  In 
order  to  define  the  distorted  body  geometry,  variables  are  introduced 
which  define  the  centerline  displacement  R,,,  the  meridional  direction  of 
the  displacement  and  the  angular  rotation  of  the  radius  vector  from  its 
original  position  For  simplicity,  the  restriction  is  made  that  the 

warpage  of  the  body  centerline  lies  in  one  plane;  therefore,  the  meridional 
direction  of  the  displacement  ^ is  the  s^aie  alonfc  the  length  of  the  body. 


SINGLE  BODY  SOLUTION 

The  representation  of  the  flow  about  a single  isolated  body  at  angle 
of  attack  is  accomplished  by  superimposing  the  axisymmetric  solution  and  the 
cross-flow  solution.  Details  of  this  analysis  are  discussed  in  the  follow- 
ing sections. 


Axisymmetric  Considerations 


Since  the  axis  of  the  body  is  not  necessarily  aligned  with  the  flow, 
only  the  component  of  freestream  velocity  in  the  axial  direction  of  the  body 
is  used  for  the  axisymmetric  solution.  The  solution  technique,  as  given  in 
detail  in  Reference  2, consists  of  distributing  a number  of  point  sources 
along  the  body  axis  at  equally  spaced  intervals  and  finding  their  strengths 
hy  simultaneous  solutions  such  that  the  body  boundary  condition  relating  the 
trelocity  with  the  slope,  S,  of  the  body  surface  is  satisfied.  The  mathemat- 
ical expression  of  this  boundary  condition  is 
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U + U 
x x. 


- S 


(3) 
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IVXTV'*  *«-“  coefficients  f„r 

ources  can  be  expressed  as  ’ undary  condition  equation  for  Ns 
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i=i 


{(C  - c 

r x 


s)  m } = s U 

1 x 


(4) 


fiedC°ntr01  P°intS  116  is  irrevalen^nd  V ^ meridi°nal  P^e  in  which 
ied  all  -round  the  body,  it  has  k f tne  boundary  condition  is  satis 

Points  should  be  take-i  L *1/  haS  b en  found,  however  that  atls“ 

hasr  cake.i  at  the  same  axial  ct,  t ’ tflat  th«  control 

greement  of  experiment  and  theory.  011  38  the  point  sources  for 


Cross-Flow  Considerations 

cross-flow  sol  lit-  -f  r\r*  j 

sS-r1™:  x/sss- -a?*  -s- 

(Uv  ) ls\Tr*triC  °r  Warped  body  cases  ther°SriSeCtl0nal  plane>  • For 
body  axis  rf  38  the  freesbream  velocity  comnon^  Vel°Clty  component 
(i.e,  normal  to  the  cross-sectional  plane)'  ParaUel  t0  the  local 

^Ith  the  single  body^ource'at^hat ' 'st t"  ltions  coincident  in  location 

ii  uycr:0srsectionai  pia-  -hia:  o{  the  ~ 

r :b 

s ssr~ 

ce.  This  can  be  expressed  as 


Ur  + Ur  + U 
® r_  r 


- S 


(5, 


Rearranging  Equation  (5)  yields 
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(6) 


5 


which  can  be  rewritten  in  terms  of  velocity  influence  coefficients  to 
yield 


N 


I {(C  - C S)  y.}  = (U  + U ) S - (U  + U ) (7) 

1 . X • 1 X X 1 I 

d d <o  s to  g 

i=l  1 

Applying  Equation  (7)  to  Ns  control  points  on  the  body  yields  the  necessary 
number  of  simultaneous  equations  for  the  solution  of  the  doublet  strengths. 
To  avoid  exaggeration  of  the  doublet  strengths,  the  control  points  for 
their  solution  are  located  such  that  the  effect  of  the  doublets  on  the 
radial  velocity  at  the  control  points  is  maximized.  Therefore,  the  control 
points  are  located  at  the  intersections  of  the  doublet  axes  and  the  body 
surface. 


TWO  BODY  SOLUTION 

The  purpose  of  the  two  body  solution  is  to  find  a set  of  source-sink 
pairs  and  a set  of  doublets  in  each  of  the  two  bodies  that  will  preserve 
the  body  shapes  when  the  two  bodies  are  in  close  proximity.  These  are 
called  the  image  source-sink  pairs  and  image  doublets,  respectively.  The 
first  step  is  the  solution  of  the  image  source-sink  pairs  while  ignoring 
the  effects  of  the  single  body  doublets  and  freestreara  cross-flow.  This 
solution  is  a modification  of  the  two  body  solution  for  two  similar  aligned 
bodies  which  was  done  in  the  work  described  in  Reference  2.  The  two- 
dimensional  Milne-Thomson  (Reference  4)  circle  theorem  (see  Appendix  B)  is 
used  to  determine  the  geometry  for  the  image  system  of  three-dimensional 
source-sink  pairs  that  are  located  at  each  station  of  the  bodies. 


Source-Sink  Image  Solution 

As  indicated  by  the  circle  theorem,  the  sink  of  the  source-sink  pair 
at  each  station  is  located  on  the  axis  (coinciding  with  the  location  of  the 
single  body  source  at  that  station)  while  the  source  is  displaced  a distance, 
6,  perpendicular  to  the  axis  in  the  direction  of  the  other  body.  The 
magnitude  of  6,  as  derived  in  Appendix  B,  is  given  by  the  equation 


The  strengths  of  the  2NS  image-pairs  (Ns  in  each  body)  are  found  by 
the  solution  of  2NS  simultaneous  equations  using  the  tangent  flow  boundary 
condition  which  is  applied  at  a control  point  on  the  body  surface  at  each 
station  of  the  two  bodies.  The  control  point  at  each  station  is  selected 
such  that  it  is  colinear  with  the  source  and  sink  of  the  image  pair.  The 
mathematical  expression  for  the  boundary  condition  is 


U + U 
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Rearranging  Equation  (9)  yields 


U-U  S - (U  + U ) 
r x x x 

p p 00  s 


(9) 


(10) 


Using  the  velocity  influence  coefficients  as  derived  in  Appendix  A, 
Equation  (10)  becomes 
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Rearranging  Equation  (11)  yields 
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Applying  Equation  (12)  to  the  control  points  at  each  station  of  both  bodies 

yields  the  required  number  of  equations  (2N  ) discussed  previously. 
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Doublet  Image  Solution 

An  image  system  must  be  derived  that  will  preserve  the  body  shapes 
when  the  effects  of  the  doublets  of  both  bodies  are  included.  Again 
the  two-dimensional  approach  is  taken.  The  image  system  that  creates  a 
circular  streamline  in  the  flow  field  of  a two-dimensional  doublet  is 


a 


derived  in  Section  8.81  of  Reference  4.  The  system  consists  of  a doublet 
located  a distance  a^/f  from  the  center  of  the  circle  on  the  line  connect- 
ing the  center  of  the  circle  and  the  original  doublet  (See  Figure  3).  The 
orientation  of  the  image  doublet  is  such  that  its  axis  is  antiparallel  to 
the  axis  of  the  original  doublet.  The  strength  of  the  image  doublet,  y', 
is  given  by 

2 

(13) 

If  an  additional  circle  that  is  centered  about  the  original  doublet 
is  considered,  then  an  image  doublet  is  needed  therein  to  counteract  the 
effects  of  the  image  doublet  in  the  first  circle  (the  effects  of  the 
original  doublet  on  the  second  circle  do  not  need  to  be  considered  since 
they  are  necessary  to  counteract  the  freestreaa  cross-flow).  Conversely, 
another  image  doublet  is  needed  in  the  first  circle  due  to  the  effects  of 
the  image  in  the  second,  and  so  on  until  a limiting  condition  is  met  (See 
Figure  4) . 

The  magnitudes  of  the  image  doublet  displacements  follow  the 
same  series  of  mathematical  expressions  as  that  of  the  source  of  the  image 
source-sink  pairs,  with  the  limiting  displacement  for  the  infinite  iteration 
given  by  Equation  (8).  Thus,  the  image  system  for  two  circles  with  one 
centered  about  a doublet  consists  of  an  infinite  number  of  image  doublets 
in  both  circles.  The  doublets,  in  the  circle  that  does  not  enclose  the 
original  doublet,  will  all  have  the  same  orientation  (antiparallel  to  the 
original  doublet)  and  will  be  displaced  from  the  circle  center  by  amounts 
ranging  from  the  initial  image  location  a2/f  to  the  limiting  value  6a 
given  as 


t 


b^  - a2 


6 - — 

a 


b2  - a2 


The  image  doublets,  in  the  circle  that  is  centered  about  the  origin? 
doublet,  will  all  have  an  orientation  that  is  parallel  to  the  original 
doublet  and  will  be  displaced  by  amounts  ranging  from  b2/(f-a2/f)  as 
shown  in  Figure  4 to  the  limiting  value  5.  given  as 

D 


a2-b2 


- 


If  an  additional  doublet  is  considered  that  is  located  in  the  first  circle, 
then  an  additional  set  of  an  infinite  number  of  image  doublets  are  needed 
in  both  circles  (See  Figure  5) . 


For  simplification,  one  image  doublet  is  used  in  each  cross-section 
to  approximate  the  infinite  number  of  doublets  suggested  by  the  two- 
dimensional  analysis.  This  cross-sectional  location  of  the  doublet  was 
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chosen  by  computational  experimentation  as  the  midpoint  between  the  first 
and  third  order  terms;  that  is,  between  a^/f  and  a2/(f-b2/(f-a2/f )) . The 
orientation  of  the  image  doublet  is  given  by  the  vector  sum  of  the  first 
and  second  order  image  terms  in  that  cross-section.  This  is  justified 
by  use  of  the  complex  potential  for  doublets  in  two-dimensional  flow.  The 
general  complex  potential  for  a doublet  as  derived  in  Section  8.23  of 
Reference  4 is 

in 


ye 


z-z 


(16) 


Therefore,  if  a single  doublet  is  used  to  represent  two  doublets  that  are 
located  at  the  same  point  but  have  different  strength  and  orientations, 
then  the  relationship  between  the  complex  potentials  must  be: 


Z-Z 


in. 

yje  1 

Z-Z 


in0 
y2e  2 

z-z 


(17) 


Multiplying  through  by  Z-Z  and  writing  the  results  in  the  form  of  cosines 
and  sines,  gives 

y^(cos  n^  + i sin  n^) 


u^(cos  n^  + i sin  n^  + y2(cos  n2  + i sin  n2) 


(18) 


By  setting  the  real  and  imaginary  parts  equal.  Equation  (18)  yields  the 
following  simultaneous  equations 


y^  cos  • Wj  cos  n^  + u2  cos  n2 


sin  h-j  ■ • i i sin  n^  + y2  sin  n2 


(19a) 

(19b) 


Dividing  Equation  (19a)  into  (19b)  yields 


tan  n,  ■ 

3 y 


y^  sin  n^  + u2  sin  n2 


^ cos  + y2  cos  n2 


(20) 


Therefore,  is  the  orientation  of  the  results  ol  the  vectorial 
addition  of  the  complex  vectors  y^eir,l  and  y2e*n2. 

Applying  Equation  (20)  to  the  combination  of  the  first  and  second 
image  doublets  yields 
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tan  n,  * 


{pj*  sin(n-n,  ) + p"  sin  n } 
b b a a_ 

{p'  cos(n-ri,  ) + p"  cos  n } 
b baa 


(21) 


The  strengths  of  the  first  and  reflected  image  doublets,  p^  and  p"  , 
respectively,  are  given  by  the  two-dimensional  analysis  as 


and 


(22a) 


(22b) 


The  two-dimensional  analysis  for  the  image  doublets  is  applied  to  the 
three-dimensional  flow  problem  in  the  same  fashion  as  the  application  of 
the  two-dimensional  analysis  for  the  sources,  where  Ra,  Rt  and  Y0  cor- 
respond to  a,  b,  and  f,  respectively.  By  these  means,  there  is  an  .'.mage 
doublet  located  in  the  cross-section  at  each  of  the  Ns  stations  of  both 
bodies . 


The  strengths  of  the  image  doublets  are  found  by  using  the  tangent 
flow  boundary  condition  which  is  now  expressed  as 


U + U + U + U + U 
r r r r . r, 

°°  s p d 0 

u +u  + u + u + u 

x x x x.  xt 

°°  s p d 6 


Rearranging  Equation  (23)  yields 


(23) 


(24) 


Expressing  Equation  (24)  in  terms  of  influence  coefficients  of  the 
doublets,  gives 


2 N 

s 


Z Z 

k-1  i-1 


- C S)p 

ik  \k 


(U  + 


u + 

x 

s 


U + 
x 

P 


U )S 


- (U  +u  ♦ u + u ) 

• s p rd 


(25) 
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To  avoid  exaggeration  of  the  doublet  strengths,  it  is  desirable  to 
locate  the  control  points  for  their  solution  such  that  the  effect  of  the 
doublets  on  the  radial  velocity  at  the  control  points  is  maximized.  This 
dictates  that  a control  point  is  located  at  the  cross-section  in  which  each 
doublet  lies.  Furthermore,  it  is  desired  that  the  meridional  location  of 
the  control  point  for  a given  cross-section  should  maximize  the  radial  ve- 
locity induced  at  the  control  point  by  the  doublet  in  that  cross-section. 
However,  the  calculations  required  to  find  this  location  are  excessive  and 
to  save  computer  time  an  alternate  location  is  used.  The  location  used  ia 
the  closest  (to  the  doublet)  of  the  two  intersections  of  the  doublet  axis 
and  the  body  surface.  This  control  point  is  shown  for  a representative 
cross-section  in  Figure  6.  This  location  insures  that  the  doublet  will 
have  a strong  effect  at  the  control  point  even  though  the  effect  is  not 
maximized. 


MULTIPLE  BODY  SOLUTION 

The  representation  of  a multiple  body  interference  flow  field  is 
acconplished  by  executing  the  single  body  solution  for  each  body,  then 
executing  the  two  body-solution  for  each  body  pair  in  turn,  and  finally, 
combining  the  results.  As  an  illustration  of  this  combination,  a repre- 
sentative cross-section  of  three  bodies  is  shown  in  Figure  7.  Points  a, 
b,  and  c represent  the  single  body  doublets  of  the  three  bodies.  The 
orientation  of  the  doublets  and  the  single  body  sources  that  are  also  a 
result  of  the  single  body  solution  are  not  shown  for  clarity.  Points  ab 
and  ba  represent  the  locations  of  the  image  doublets  that  are  created  by 
the  execution  of  the  two-body  solution  for  Bodies  A and  B.  Again,  for 
clarity,  the  doublet  axes  and  the  image  source-sink  pairs  that  are  also 
a result  of  the  two-body  solution  are  not  shown.  Similarly,  points  ac 
and  ca  represent  the  locations  of  the  images  created  by  the  two-body  solu- 
tion of  Bodies  A and  C and  points  be  and  cb  represent  the  locations  of 
the  images  created  by  the  two-body  solution  of  Bodies  B and  C. 

By  combining  the  results  of  all  the  two-body  solutions,  there  are 
image  systems  in  each  body  that  prevent  that  body's  shape  from  being 
distorted  by  the  single  body  sources  and  doublets  of  the  other  bodies. 

This  does  not,  however,  prevent  the  body  shape  from  being  distorted  by 
the  image  systems  of  the  other  bodies.  To  illustrate  this.  Body  A in 
Figure  7 is  used  as  an  example.  The  image  systems  created  by  the  two- 
body  solution  of  Bodies  A and  B (partially  represented  by  points  ab  and 
ba)  prevent  the  distortion  of  Body  A by  the  single  body  sources  and 
doublets  of  Body  B and  also  prevent  the  distortion  of  Body  B by  the 
single  body  sources  and  doublets  of  Body  A.  Similarly,  the  image  systems 
created  by  the  two-body  solution  of  Bodies  A and  C (represented  by  points 
ac  and  ca)  prevent  mutual  distortion  of  the  shapes  of  Bodies  A and  C. 

The  image  systems  presented  thus  far  are  required  in  order  to  prevent 
the  distortion  of  the  chape  of  Body  A by  the  single  body  sources  and 
doublets  of  Bodies  B and  C.  However,  the  image  systems  that  are  created 
by  the  two-body  solution  of  Bodies  B and  C are  also  included  in  the  flow 
field.  Hence,  an  additional  image  system  is  needed  in  Body  A (partially 
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Doublet  Axis 


Figure  6.  Example  of  Control  Point*  Used  for 
Solution  of  Image  Doublet  Strengths 
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represented  by  points  abc  and  acb)  to  prevent  distortion  by  the  image 
source-sink  pairs  and  image  doublets  resulting  from  the  two-body  solution 
of  Bodies  B and  C.  A similar  situation  exists  for  Body  B and  for  Body  C 
due  to  the  two-body  solutions  of  Bodies  A and  C and  Bodies  A and  B, 
respectively.  It  has  been  found,  however,  that  the  contribution  of  these 
higher  order  image  terms  is  at  best  an  order  of  magnitude  less  significant 
than  the  two-body  image  system.  Therefore,  since  these  terms  have  little 
effect  on  the  solution,  as  can  be  shown  by  the  fact  that  the  distortion 
of  the  body  shapes  is  small  when  only  the  results  of  the  single  body  and 
two-body  solutions  are  added,  the  higher  order  image  terras  are  neglected. 
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SECTION  III 


NUMERICAL  EXAMPLES 


The  numerical  examples  presented  here  are  based  on  the  M-117  bomb 
shape  (without  fins)  which  is  shown  in  Figure  8. (a) . The  geometry  of  the 
M-117  bomb  can  be  described  in  three  axial  intervals.  They  are  the  nose 
section,  the  center  section,  and  the  tail  section.  The  nose  section  is  of 
ogival  shape  (circular  arc)  with  a length  to  maximum  diameter  ratio  of 
1.32.  Tha  center  section  is  a circular  cylinder  with  a length  of  1.22 
diameters.  The  tail  section  is  a straight  taper  of  7.5  degrees  and  termin- 
ates in  a boatt-ail.  The  length  of  the  actual  M-117  tail  is  2.96  diameters; 
however,  a length  of  3.47  diameters  is  used  in  these  examples  to  give  some 
ullowance  for  flow  separation  at  the  truncated  boattail. 

All  theoretical  data  presented  in  this  section  were  obtained  by  the 
use  of  a Fortran  Computer  Program.  The  pressure  coefficients  along  the 
body  surface  were  calculated  for  24  evenly  spaced  meridional  planes. 

Force  and  moment  coefficients  were  then  calculated  by  numerically  integrat- 
ing the  pressure  coefficients.  The  reference  area  and  length  used  for  the 
coefficients  are,  respectively,  the  maximum  cross-sectional  area  and  the 
diameter  of  the  M-117  bomb.  The  moment  coefficients  are  referred  to  the 
center  of  gravity  which  is  2.05  bomb  diameters  from  the  nose.  To  aid  in 
the  evaluation  of  the  modeling  technique,  error  calculations  were  made  at 
each  station  around  the  body  by  finding  the  flow  angle  (with  respect 
to  the  body  centerline)  created  by  all  induced  and  freestream  velocities 
and  subtracting  the  angle  of  the  body  surface. 

The  configurations  used  for  the  numerical  examples  are  shown  in  Figvre 
9.  They  are  drawn  as  though  the  viewer  were  looking  downstream  at  the 
nose  of  the  bodies.  The  arrows  indicate  the  direction  of  the  freestream 
cross-flow  due  to  angle  of  attack. 


SINGLE  BODY 

The  case  of  a single  M-117  bomb  aligned  with  the  flow  has  been  pre- 
viously treated  in  Reference  2, and  the  resulting  pressure  distribution  has 
been  found  to  be  in  good  agreement  with  experimental  data  obtained  in  the 
low  speed  wind  tunnel  at  the  contractor’s  as  shown  in  Figure  10.  Figure 
11  shows  the  pressure  distributions  along  the  leeward  and  windward  meridi- 
onal planes  for  a M-117  bomb  at  5 degrees  angle  of  attack. 

The  analytical  force  and  moment  coefficients  for  the  isolated  body 
were  found  to  be  negligible  except  for  the  pitching  moment  which  is  plotted 
in  Figure  12  for  angles  of  attack  ranging  from  zero  to  10  degrees. 
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Figure  8.  Profiles  of  M-117  Bomb  Shapes 


|Cross-Flov 


Cross-Flow 


(a)  Single  Body 


) Two  Body 


NOTE:  All  dimensions  in 

c'  ^r«e  Body  body  diameters 


Figure  9. 


Sin“'  ,n<i  »«dy  Configurations 


(Carriage  Position) 
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ion  - Single  Body 


Error  calculations  for  a single  body  at  angle  of  attack  showed  that 
the  difference  between  the  inclination  of  the  streamlines  at  the  body 
surface  for  points  in  meridian  planes  other  than  that  containing  the  con- 
trol points  and  the  inclination  of  the  actual  body  surface  did  not  exceed 
0.0005  degree  at  any  of  the  points  checked. 

TWO  PARALLEL  BODIES 

The  case  of  two  parallel  M-117  bombs  aligned  with  the  flow  was  treated 
in  Reference  2.  As  in  Reference  2,  the  distance  between  the  center  lines  of 
two  bodies  used  for  this  case  is  1.042  diameters.  This  corresponds  to  the 
separation  of  two  M-117  bombs  when  mounted  on  a Triple  Ejector  Rack  (TER). 

The  resulting  pressure  distribution  was  found  to  be  in  good  agreement  with 
experiment  as  shown  in  Figures  13  and  14.  No  experimental  data  are  avail- 
able for  this  two-body  case  at  an  angle  of  attack;  however,  the  relative 
change  in  the  pressure  distributions  is  shown  in  Figure  15.  As  can  be  seen 
from  this  figure,  the  pressure  coefficient  at  the  point  between  the  bodies, 

9=180  degrees,  is  apparently  insensitive  to  the  angle  of  attack.  However, 
there  does  seem  to  be  a significant  change  with  angle  of  attack  at  0=0  degrees. 

The  effect  of  angle  of  attack  on  side  force  and  pitching  moment  coefficients 
is  shown  in  Figures  16  and  17. 

§ 

To  explore  the  validity  of  this  image  system  a check  was  made  on  the 
body  boundary  conditions  at  meridian  planes  other  than  that  containing  the 
control  points.  The  poorest  match  of  the  tangent  flow  boundary  condition 
occurred  along  the  outside  meridian  and  the  magnitude  of  the  error  incurred 
was  unaffected  by  the  angle  of  attack.  A plot  of  the  difference  between 
the  body  surface  angle  and  the  flow  angle  along  this  meridian  is  shown  in 
Figure  18. 

THREE  B0DIZS 

The  three  body  case  represents  three  M-117  bombs  in  their  relative 
positions  as  if  mounted  on  a TER  and  is  shown  in  Figure  9(c).  Experimental 
data  for  this  configuration  at  zero  angle  of  attack  were  obtained  from  a 
pressure  model  and  two  models  of  the  M-ll^  bomb  in  the  low-speed  wind 
tunnel  of  the  contractor.  The  results  of  the  wind  tunnel  test  and  theory 
for  the  top  meridional  plane  of  the  central  body  are  shown  in  Figure  19. 

As  can  be  seen,  the  agreement  between  theory  and  experiment  is  poor  along 
the  aft  portion  of  the  body.  Subsequent  wind  tunnel  tuft  tests  have 
shown  that  there  is  flow  separation  between  the  bodies  along  the  tails  of 
all  three  bodies.  It  was  for  this  reason  that  it  was  decided  to  analytically 
modify  the  shape  of  the  M-117  to  try  to  account  for  flow  separation  along 
one  side  on  the  tail  of  the  body. 

The  non-axisyranetric  configuration  consideied  consists  of  the  normal 
M-117  nose  and  center  sections.  The  tail  section,  however,  has  a decreased 
taper  and  angled  centerline  such  that  one  side  maintains  the  original 
contour  of  the  M-117  tail  while  the  opposite  side  is  displaced  [See  Figure 
8(b)].  All  three  bodies  were  modified  in  this  manner  with  the  tail  of  the 
central  body  being  warped  upward  and  the  tails  of  the  two  upper  bodies 
being  warped  toward  the  central  body.  The  center  line  of  the  tail  of  the 
modified  body  was  set  at  an  angle  of  2 degrees.  The  pressure  distribution 
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Figure  14.  Pressure  Distribution  - Two  Bodies  at  0-Degree  Angle  of  Attack  (0=0°) 


Comparison  of  the  Analytical  Pressure  Distributions  for  Two  Bodies 
at  0-  and  5-Degree  Angle  of  Attack 
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Degree  Angle  of  Attack  (Concluded) 


PITCHING  MOMENT  COEFFICIENT 


figure  17.  Pitching  Moment  Coefficient 


Two  Body  Configuration 
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Figur*  IS.  Boundary  Condition 


SECTION  IV 


CONCLUSIONS 


The  force  and  moment  coefficients  plotted  in  the  numerical  examples 
show  significant  variation  with  angle  of  attack;  linear  in  the  case  of  two- 
body  pitching  moment,  nonlinear  in  the  case  of  two-body  side  force.  This 
indicates  that  for  complete  investigation  of  interference  effects  the 
freestream  cross-flow  due  to  angle  of  attack  should  be  taken  into  account. 

The  modification  of  body  geometry  to  account  for  flow  separation  was 
found  to  be  Inadequate  for  the  particular  case  tried  in  the  numerical 
examples.  It  was  found,  however,  that  some  of  the  characteristics  of  the 
separated  flow  were  simulated  by  the  distorted  body  and  it  is  speculated 
that  this  method  might  be  used  successfully  for  flow  separation  problems 
that  are  less  complex  than  the  one  tried  in  the  numerical  examples. 

It  is  also  suggested  that  the  capability  of  modeling  slightly  non- 
axlsymmetrlc  bodies  will  be  useful  for  many  flow  problems.  One  foreseen 
application  of  this  technique  is  the  representation  of  the  flow  field  in 
the  vicinity  of  an  aircraft  fuselage  by  using  multiple  non- axisymme trie 
bodies  that  approximate  the  shape  of  the  fuselage. 
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APPENDIX  A 


INFLUENCE  COEFFICIENTS 


Source  Flow 

The  potential  function  of  a three-dimensional  source  as  given  in  Section 
15-20  of  Reference  4 is  $ * m/r  (reversing  the  sign  convention  results  in 
$ “ - m,/r).  This  function  can  be  described  in  a three-dimensional  coordinate 
system  with  the  source  located  at  the  origin  by  letting  r = (x2+y2+z2)*S; 
therefore 


♦ - T (A-l) 

(x2+y2+z2)*5 

The  velocities  in  the  x,  y,  and  z directions  that  are  induced  by  the 
source  at  the  point  (x,y,z)  can  be  found  by  the  partial  differentiation  of 
the  field  potential: 
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The  velocity  influence  coefficients, 
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and  C are  defined  as 
z 

s 


(A-5) 


and 


m 


(A-6) 


(A- 7) 


37 


i . ■ iifwnwim  i 


Therefore,  using  Equations  (A-2,  A-3,  and  A-4),  then 


Xs  {x2+y2+z2} 3^2 


(A-8) 
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S {x2+y2+z2  }3/2 


(A-10) 


Doublet  Flow 

The  potential  function  of  a three-dimensional  doublet  of  strength  y is 
given  in  Section  15.26  of  Reference  4 (after  sign  reversal  again)  as 

cos  0 

♦ " - V 7T 

If  the  doublet  is  considered  at  the  origin  of  an  x,  y,  z axis  system  with 
the  doublet  axis  oriented  in  the  negative  y direction,  then  r becomes 
(x2+y2+z2)1/2,  cos  0 becomes  - y/(x*+y2+z2)l/2  and  the  velocity  potential 
becomes 


,2  -,2.. 2n 3/2 


(A-ll) 
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Differentiating  Equation  (A-lly  yields  the  velocity  components: 
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The  velocity  influence  coefficients 


U 

z 


of  the  doublet  are  defined  as 

(A-15) 

(A-16) 

(A-17) 


Using  Equations  (A-12,  A-13,  and  A-14),  then 
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Transformation 

The  velocity  coefficients,  of  course,  can  be  transformed  into  another 
coordinate  system  if  desired.  For  example,  if  a source  of  strength  m is 
located  inside  of  a circle  as  shown  in  Figure  A-l,  with  the  x axis  point- 
ing into  the  paper,  then  the  radial  and  tangential  velocities  (U,  and  Ufc) 
induced  by  the  source  m,  at  the  point,  P,  can  be  found  by  transforming 
the  velocity  components  in  the  y and  z directions  through  the  angle  B. 

7b.!  3 gives 

U - U cos  B + U sin  B (A-21) 

r y z 

Ufi  - U cos  B - U sin  B (A-22) 

o z y 


U^ing  the  velocity  influence  coefficients  as  previously  derived,  then, 

U ■ C a cos  B + C m sin  B (A-23) 

r y z 

8 e 

U ■ C m cos  B - C m sin  B (A-24) 

t z y 

8 's 


40 


Defining  the  radial  and  tangential  velocity  influence  coefficients  of  the 
source  as 


and  C = — , then 

t m ’ 
s 

Cr  “ C cos  0 + C sin  B (A-25) 

S yB  3 

and 

Ct  " Cz  C0S  6 " Cv  8in  6 (A"26> 

s s ys 

Comparing  Equations  (A-25  and  A-26)  with  Equations  (A-21  and  A-22) , 
it  can  be  seen  that  the  velocity  influence  coefficients  can  be  transformed 
as  if  they  were  actual  velocities. 


APPENDIX  E 


IMAGE  SYSTEM  FOR  TWO  CIRCULAR  CROSS-SECTIONS 
IN  TWO-DIMENSIONAL  FLOW 


The  problem  discussed  here  has  been  previously  treated  in  Section  IV 
of  Reference  1 and  is  discussed  here  as  an  aid  to  the  reader. 

From  the  single  body  solution  for  an  axisymmetric  body,  there  exist  a 
number  of  point  sources  located  along  the  body's  axis  that  are  vised  to 
satisfy  the  tangent  flow  boundary  condition  at  the  body  surface.  If  two 
such  mathematical  models  are  brought  into  close  proximity,  then  the  sources 
in  each  body  will  create  perturbations  at  the  other  body's  surface  and 
therefore  will  be  mutual  distortion  of  the  body  shapes.  A two-dimensional 
approach  is  taken  to  find  an  appropriate  image  system  that  will  preserve 
the  body  shapes. 

A cross-section  of  the  two  bodies  is  represented  in  two-dimensional 
flow  by  two  circles  with  a source  at  each  center  as  shown  in  Figure  B-l. 

It  is  desired  that  the  circular  shape  of  Body  A be  maintained  by  creating 
an  image  system  that  will  cancel  the  perturbations  induced  by  the  source 
in  Body  B that  are  normal  to  the  surface  of  Body  A.  The  image  (derived 
in  Section  8.61  of  Reference  4)  for  a source  outside  of  a circular  cylinder 
satisfies  this  condition.  Therefore,  the  image  system  required  consists  of 
a sink  of  strength  m at  the  center  of  A and  a source  of  equal  strength 
displaced  from  the  center  of  A toward  the  center  of  B by  an  amount  a2/f 
where  a is  the  radius  of  A and  f is  the  distance  between  the  centers  of  A 
and  B. 

An  image  system  is  also  desired  in  B to  counteract  the  source  of  A. 

This  will  consist  of  a sink  of  strength  ma  at  the  center  of  B and  a source 
of  equal  strength  displaced  from  the  center  of  B toward  the  center  of  A 
by  an  amount  bvf  where  b is  the  radius  of  B (See  Figure  B-2). 

Additional  images  are  now  needed  in  A and  B to  counteract  the  effects 
of  the  additional  sources  and  sinks  added  to  A and  B.  That  is,  a sink  of 
strength  ma  at  the  center  of  A and  a source  of  equal  strength  displaced 
a distance  a?/(f-b2/f)  from  the  center  of  A are  needed  to  counteract  the 
effect  of  the  image  source  in  B.  Also,  a source  of  strength  ma  at  the 
center  of  A and  a sink  of  equal  magnitude  displaced  from  the  center  of  A 
by  an  amount  a2/f  are  needed  to  counteract  the  effects  of  image  sink  in  B. 

A similar  situation  is  found  in  B.  After  these  images  are  added,  then 
additional  images  are  needed  for  them  also,  and  so  on  to  infinity.  Figure 
B-3  shows  the  results  of  several  iterations. 
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The  displacement  of  the  sources  from  the  center  of  the  circle  reaches 
a limiting  point  (See  Figure  B-3)  and  the  displacement  of  the  infinite 
iteration  is  expressed  as  (using  the  displacement  in  circle  A as  ar  examp  e) 

6 = a2/(f-b2/(f-a2/(f-b2/(f-a2/(. .,.)))) 

(B-l) 


Examining  Equation  (B-l),  it  can  be  seer,  that  due  to  the  infinite  series 
the  equation  can  be  rewritten  as 

ar 

6 = f“  b^~  (B— 2) 

f-6 


Rearranging  Equation  (B— 2) , it  can  be  seen  that 
(b2-a2^ 


f _ (b2-a2) 
" _ f 


6 + az  = 0 


(B-3) 


Solving  the  quadratic  equation  for  6, 

1 [ b2-a2l 

* ■ i Lf  - — J 


(B-4) 


(The  solution  with  a positive  radical  lies  outside  the  circle  and  is 
therefore  not  considered) . 


A similar  solution  can  be  found  for  the  limit  of  the  displacement 
in  B and  the  result  is  the  same  as  Equation  (B-4)  except  that  a and  b are 

switched. 
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LIST  OF  SYMBOLS 


a Local  radius  of  Body  A 

A Denotes  Circle  A or  Body  A 

b Local  radius  of  Body  B 

B Denotes  Circle  B or  Body  B 

C Pitching  moment  coefficient 
m 

C Pressure  coefficient 
P 

C^  Radial  velocity  influence  coefficient  (at  a control  point) 
C Side  force  coefficient 

C^  Axial  velocity  influence  coefficient  (at  a control  point) 

f Distance  between  the  centers  of  two  circles 

m Source  strength  is  strength  of  ith  source) 

N Number  of  stations  along  body 
s 

R Radius  of  body  cross-section 
R^  Displacement  of  body  centerline 
S Slope  of  body  surface  with  respect  to  centerline 
Radial  velocity 
Axial  velocity 

V Velocity  vector 

w Complex  potential  function 

x,y,z  Rectangular  coordinates  (x  along  body) 

Y Distance  between  centerline  of  Body  A and  Body  B measured 
perpendicularly  to  centerline  of  Body  A 
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LIST  OF  SYMBOLS  (CONTINUED) 


Z 

6 

6 


n 

e 

u 

* 

Tv 

* 

V 

V:1 


Complex  coordinates 
Transformation  Angle 
Source  displacement  distance 
Orientation  of  doublet  axis 
Angular  coordinate 
Doublet  strength 
Velocity  potential 

Angular  location  of  displacement  of  body  centerline 

Rotation  of  body  cross-section  due  to  displacement  of  body  centerline 

Laplacian  operator 

Source  location 

Doublet  Location 
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LIST  OF  SYMBOLS  (CONCLUDED) 


Subscripts 

a Denotes  a relationship  with  Body  A 

b Denotes  a relationship  with  Body  B 

d Indicates  a property  associated  with  a doublet 

i,k  Refers  to  the  i1"^  source  or  the  iC^  doublet  or  i^  body  or  k* 

body. 

p Denotes  a property  associated  with  image  pairs 
v Denotes  radial  property  or  value 
s Indicates  a property  associated  with  a source 

t Indicates  tangential  property  or  value 

x Denotes  axial  property  or  value 

® Denotes  freestream  value 

6 Denotes  property  associated  with  displaced  or  image  doublets 
Superscripts 

+ Indicates  property  associated  with  source 
(+  implied  if  omitted) 

Indicates  property  associated  with  sink 
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